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Abstract 

We use the numerical renormalization group to calculate the auxiliary spec- 
tral functions of the U = oo Anderson impurity model. The slave-boson and 
pseudo-fermion spectral functions diverge at the threshold with exponents af, 
and atf given by the X-ray photoemission and the X-ray absorption expo- 
nents respectively. In contrast to the NCA, the exact exponents obtained here 
depend on the impurity occupation number. In general, vertex corrections in 
the convolution formulae for physical Green's functions are singular at the 

threshold and may not be neglected in the Fermi liquid regime. 
PACS numbers: 71.27. +a,71.28.+d,72.20.Hr 
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Strongly correlated fermi systems have been under intense investigation over the past 
few years. This development was triggered by the discovery of high temperature supercon- 
ductors, the subsequent proposals of an electronic mechanism for superconductivity [jlj] and 
the many fascinating ideas concerning the normal state 0. Despite much theoretical effort, 
a systematic and controlled method for dealing with such systems is still lacking. At the 
heart of the problem is the influence of a strong on-site Coulomb repulsion U on a system of 
electrons living in a lattice in d = 2 dimensions. In the limit of U — > oo one has to deal with 
a projection of the dynamics into the part of the Hilbert space without any doubly occupied 
lattice sites. It is difficult to implement this projection in conventional perturbation the- 
ory. Therefore it has been proposed some time ago to effect the projection by introducing 
auxiliary operators in a larger Fock space ||. Thus for spin | fermions in a single band 
model, where there are three states per lattice site in the U — > oo limit, one introduces a 
Bose operator b\ and a pair of Fermi operators f- a (a = f , J.) acting on their vacuum to 
create the empty lattice site i and the singly occupied lattice site i with j or |. As long as 
the number of auxiliary particles is exactly one at each site, Q = b\bi + J2a fLfiv = 1> the 
above definition constitutes a faithful representation of the original model. Unfortunately, 
the model of interacting bosons and fermions defined in this way for any physical model of 
this type (e.g. the t-J model) is difficult to solve, and few reliable results have been obtained 
so far. 

In order to develop approximation schemes it is therefore of interest to apply the auxiliary 
particle representation to simpler models, for which comparison with exact results is possible. 
One such model is the infinite U single impurity Anderson model of a magnetic moment in 
a metal |4j]. 

The physics of this model, in particular the Kondo effect, has been clarified long ago 
using qualitative scaling arguments ||, a numerically exact renormalization group (NRG) 
analysis ||, and the exact solution based on the Bethe ansatz method 0j. A quantitative 
analytic theory of the dynamical properties is not completely available yet, although the 
numerical renormalization group has been shown to yield close to exact results for both zero 
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and finite temperature dynamics |J. Despite this it is desirable to have an alternative way 
of reliably calculating dynamic properties of this model which is more easily generalized 
to lattice models. The slave-boson approach offers such a possibility. It has been studied 
within two different approximations. The first one is a mean-field theory of the bosons and 
the constraint field 0. It appears to qualitatively describe the Fermi liquid behaviour at low 
temperatures known to exist in this model. However, it can not describe the smooth cross- 
over to the free magnetic moment behaviour at high temperatures. The second one involves 
a self-consistent summation of perturbation theory and exact projection on the physical 
subspace known as the "non-crossing approximation" or NCA [IIJ. It is found to describe 
correctly the cross-over from the high-temperature to the low-temperature phase, except for 
the appearance of spurious non-analytic behaviour in the physical Green's functions below 



a characteristic temperature T NC a, where T NC a << T K , the Kondo temperature pl| , |r2 
Within the NCA the auxiliary spectral functions for pseudo-fermions (/) and slave-bosons 
(b), At b (u), above the threshold, E , behave as A^ b (uj) ~ \u — E \~ a f> b , with aj = 1/(N + 1) 
and «b = N/(N + 1), where N denotes the spin degeneracy. We note in passing that the 
corresponding NCA exponents for the M-channel Anderson model have recently been shown 
to give the correct singular behaviour [TB[] expected for M > 1. 



In order to correctly describe the low-frequency behaviour in such models, a clear un- 
derstanding of the infrared behaviour of the auxiliary spectral functions is essential. In 
this paper we address the question of the true behaviour of the auxiliary spectral func- 
tions at zero temperature and in the limit of frequency uo — > E and clarify the origin of 
the non-analyticities in the physical Green's functions within the NCA. We adapt Wilson's 
renormalization group method to the auxiliary particle problem and attempt to calculate 
the spectral functions numerically. 

The Anderson model in auxiliary particle representation takes the form 

H = H c + e d J2 flU + V + h.c), (1) 

(7 (7 

where H c = J2ka e kcl a Cka is the conduction electron kinetic energy and Co a = J2k c ka anni- 



hilates a conduction electron with spin a at the impurity site 0. Following Wilson || we 
(i) linearize the spectrum about the Fermi energy — > k, (ii) introduce a logarithmic mesh 
of k points k n = A~ n and (iii) perform a unitary transformation of the Ck a such that co a is 
the first operator in the new basis and H c takes the form of a tight-binding Hamiltonian in 

k-space, 

oo 

He = £ E^ A ~ n/2 ( C n+l<x C - + ^ C -), ( 2 ) 
n=0 

with £ n — > (1 + A~ 1 )/2 for n >> 1. These steps are explained in detail in || and can be 
taken over for the present model without change. 

The Hamiltonian ([!]) together with the discretized form of the kinetic energy (0) in the 
new basis is now diagonalized by the following iterative process: (i) one defines a sequence 
of finite size Hamiltonians H N by replacing H c in (|) by = J2n=o So- Cn^~~ n ^ 2 ( c n+ia c n<T + 
h.c); (ii) starting from H = edJ2<j ftf<r + VJ2a( c oir^fa- + h.c.), each successive hopping may 
be considered as a perturbation on the previous Hamiltonian; (iii) the Hamiltonians are 
scaled such that the energy spacing remains the same. This defines a renormalization group 
transformation H N+ i = A 1 / 2 ^^ + £ N Z) -( c !v+i - c A r o- + h.c.) — Eq,n+i, with E G<N+1 chosen so 
that the ground state energy of H^ + x is zero. The diagonalization of the Hamiltonians 
is greatly simplified by the existence of the following conserved quantities: electron number 
N e , total spin S, spin projection S z and auxiliary particle number Q. The eigenstates of 
H N may therefore be labeled as \N e , S, S z , Q; k > with k denoting the remaining quantum 
numbers, and the Hamiltonian matrices have block diagonal structure. Considering that the 
smallest energy scale probed by H N is uj^ ~ A -(Af ~ 1)/2 , for a typical Kondo temperature of 
Tk which is << 10~ 2 of the conduction band width, it is necessary to iterate beyond ~ 15 
for A 2 to get below the Kondo scale. Since the dimension of is growing as 4^ it is 
necessary to truncate the higher energy states for N > 7. Wilson presents arguments why 
this is a valid approximation. We have kept the lowest 250 states in each iteration step and 
checked that our results remain unchanged when this number is increased to 2000 states. 



In Coleman's slave-boson formulation of the mixed valence problem |10| the retarded 



Green's functions for pseudo-fermions and slave-bosons are defined in the enlarged Hilbert 
space consisting of the disjoint subspaces Q — 0, 1, ... by 

GfAu, T, A) = Zah £ |M^| 2 (e-^+^ +1 » + c -/*»+*»)/ (a , - (e m - e n )), 

Q,m,n 

where —A is a chemical potential associated with the auxiliary particle number Q, Zgc{T, A) 
is the grand-canonical partition function and M^ n =< Q + 1, m\0'\Q, n > with O = 
f a , b are the many-body matrix elements for pseudo-fermions and slave-bosons respec- 
tively. We are interested in the zero temperature projected spectral functions At b = 
- lim T ^ liniA^ooflm G fjb (u, T, X)/tt] and Aj >b = - lim T ^ lin^oo e _/3w [Im G f)b (uJ, T, 

A+ b (u) = irJ2 I < 1, m|Ot|$o > \ 2 6(cu + E%f - e m ), (3) 

m 

A lM = I < $i\0% m > \ 2 6(u + e m - E%f). (4) 

m 

Here |$o > is the groundstate of the Q = subspace of non-interacting conduction electrons 
and |1, m > are the excited states of the Q = 1 subspace of the interacting system, Eqs° and 
e m are the corresponding energy eigenvalues. The spectral functions A~j b vanish identically 
below the threshold E = Eq^ 1 — Eq^ . Similarly, in the expression for Aj b , \<S>i > is the 
groundstate of the interacting system (Q = 1 subspace) and |0,m > are the excited states 
of the non-interacting conduction electron system, Eq^ 1 and e m the corresponding energy 
eigenvalues. These spectral functions vanish above the threshold energy E . 

The pseudo-particle operators /J, in (||-§D have matrix elements Mj^ n only between 
states from the subspace Q = 1 and the subspace Q = 0. These matrix elements are 
calculated for each iteration step, and are substituted together with the energy eigenvalue 
e m into (0-§D to give spectral functions A^j^u). In principle if all states up to stage 
iV were retained, would describe excitations on all energy scales from the band edge 
D = 1 down to the lowest energy scale present in H^, i.e. u^. Due to the elimination 
of higher energy states at each step, the actual range of excitations in Hn is restricted to 
lon < uj < Kojjsi, where K « 7 for A « 2 retaining 500 states per iteration. Thus at step N, 
the spectral functions are calculated at an excitation energy uo m 2u>n in the above range. 
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The delta functions in are broadened with Gaussians of width cxn ~ &n appropriate 

to the energy level structure of H^. 

The RG calculations were performed for A = 2 and keeping 250 states per iteration for 
each subspace (Q = 0, 1). The hybridization strength A = n V 2 /2D was chosen to be 0.01/} 
with the half-bandwidth D = 1. Several values of the local level position where chosen 
in order to characterize the low energy behaviour of the spectral densities in the various 
regimes. In this letter we are mainly concerned with the low frequency limit, discussion of 
high energy effects is left to a future publication. 

The auxiliary spectral functions exhibit threshold behaviour and are shown in Fig. ([!]) in 
the Kondo regime with similar behaviour near the threshold in the mixed valent and empty 
orbital regimes. We note that this is not unexpected since the projected spectral functions 
are very similar in appearance to the core-level spectral functions in the usual X-ray 



problem so one expects as a result of the orthogonality catastrophe theorem a similar 
threshold behaviour. This analogy is useful but requires care since the matrix elements in 
(Q-^) are no longer between two non-interacting systems as in the X-ray problem. This 
leads in particular to a new energy scale, T , for the onset of the asymptotic power law 
behaviour, which is Tjc, A or in the Kondo, mixed valent and empty orbital regimes 
respectively. We find that it is only in the Fermi liquid regime, \uj — Eq\ << To, that the 
power law behaviour is well characterized. 

The results for the threshold exponents are shown in Fig. (Q) as a function of n/, the local 
level occupancy at T = 0. The latter was calculated by evaluating nf(T) from the partition 
function at a sequence of decreasing temperatures T/v ~ A - ^ -1 )/ 2 and then taking the limit 
T — > 0. Remarkably these exponents turn out to be the usual photoemission and absorption 
exponents for the X-ray problem and are given in terms of the conduction electron phase- 
shift at the Fermi level, 5 a = 5 a (e F ), by 

a f = 25Jn - ]T(<Wt0 2 = n f - n)/2 (5) 

(7 

a b = l- E(^/^) 2 = 1 - n }/ 2 (6) 



where the last equations on the RHS of (0-§D follow from the Friedel sum rule, 8 a = 7r n//2. 
These results are clearly illustrated in Fig. (0) where the functions rif — nj/2 and 1 — ri 2 /2 
are plotted against n./ together with the exponents o;/^ deduced from the spectral functions. 
The exponents a/^ agree with the RHS of (|5|-§P to 3 significant figures and are the same 
below and above the threshold, 

A% = a%\u-E \- a '* (7) 

The slave-boson exponent a& corresponds to the usual XPS exponent in the X-ray problem. 
On the other hand the pseudo-fermion exponent af takes the form of an absorption exponent 
in the X-ray problem. This is at first sight unexpected since we are dealing with single- 
particle Green's functions in which the processes correspond to removing a particle from 
the system altogether. We attempt to give a qualitative explanation of this. Assuming, as 
is shown rigorously by the numerical results, that the threshold exponents are determined 
by the phase shifts we note that the removal of a slave-boson has no effect on the local 
electronic charge so the relevant phase shifts entering the XPS exponent are 8 a which are 
fixed by the Friedel sum rule. Hence the usual exponent, 1 — Y.ai^/ 71 ) 2 > appears in Af. On 
the other hand the removal of a pseudo-fermion reduces the electronic charge by one. By 
charge neutrality the sum of phase shifts Eo-^ is reduced by tt. Assuming, as in the usual 
X-ray problem Jl5| , that the total phase shift occurs in the channel a in which a fermion is 
removed we have that 1 - E^/tt) 2 -> 1 - [(6 a /n - l) 2 + (8^/n) 2 ] -> 2^/vr - £ ct (5 ct /tt) 2 
which leads to the exponent af found numerically. We note that the same functional form 
of the exponents on the phase shift (|5|-§) is also found in the spinless model with constraint. 
The single phase shift in this case is given by 8 = iruf. 

With the above results we can gain some insight into the role played by vertex corrections 
in the slave-boson method. These appear in Green's functions for slave-bosons, pseudo- 
fermions and in the convolution formulae for physical Green's functions. In general these 
are difficult to calculate within perturbative schemes but they may be estimated within 
the present non-perturbative approach. First we note that the above threshold exponents 
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generalized to the N-fold degenerate model are atf = 2n//N — nj/N and aj, = 1 — nj/N 
so the NCA exponents l/(iV+l) = 1/N + 0(1/N 2 ) and N/(N + 1) = 1 - 1/N + 0(1/N 2 ) 
are only correct in the limit rif — > 1 and iV — > oo. Away from this limit, vertex cor- 
rections in the auxiliary Green's functions, neglected in the NCA, are therefore important 
in determining the correct threshold exponents. To see the influence of vertex corrections 
on the physical Green's functions we consider specifically the impurity Green's function 
G a {ui) =« f a ; fib »q=i. This is given in terms of the projected auxiliary Green's 
functions, by 

G a (i(j) = —— ^Gf(iuo v + iuo)G b (iuo u )V(iuj u + iu, iuo u ) (8) 

P V 

where V(iu u + iu,iu v ) represents the renormalized vertex. The effect of vertex corrections 
can be seen by evaluating the impurity spectral density, p'f a (uj) = — Im [G a {uj + i5)/it}, at 
low frequencies without vertex corrections, i.e. for V — 1. Using (§-0) in (|j) gives 

p' fa {co -> 0+) = a+a^cu^f-^Bil -a f ,l + a b ) (9) 

where B is the Beta function. Thus in the absence of vertex corrections the impurity spectral 
density diverges at the Fermi level as a power law \uj\ l - a f- ab = l^l-M 1 -™/). From the Friedel 
sum rule, Pf a (uj = 0) = sin 2 (7m//2), so the vertex corrections are singular at low energies, 
i.e. close to the threshold, and lead to a singularity in pf a at the Fermi level which cancels 
that in (||). This singular behaviour together with next leading order corrections to the 
spectral functions and vertex corrections should restore the Fermi liquid behaviour which 
is well known for this model. Similar vertex corrections appear in the calculation of other 
physical quantities such as the dynamic spin susceptibility. 

To summarize, we have characterized the low energy behaviour of the auxiliary spectral 
functions of the single-channel Anderson model at T = 0. The power law behaviour close to 
the threshold is approached slowly and only sets in well below the characteristic low energy 
scale of the model. The threshold exponents were related to the conduction electron phase 
shifts and depend on the occupation number rif, in contrast to the corresponding NCA 
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exponents. Vertex corrections in the pseudo-fermion and slave-boson Green's functions are 
therefore required to obtain the correct threshold exponents. They are also needed in the 
convolution formulae for physical Green's functions to restore the well known Fermi liquid 
behaviour for the single-channel model. It is interesting to speculate whether the vertex 
corrections are non-singular at the threshold in multi-channel impurity models exhibiting 
non-Fermi liquid behaviour. Work on this is in progress and will be reported elsewhere. The 
results presented here could provide guidance in developing reliable self-consistent slave- 
boson schemes. These could then be extended to the important case of lattice models, 
which are more difficult to study within the numerical renormalization group approach. 
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FIGURES 

FIG. 1. The pseudo-fermion Aj (solid lines) and slave-boson A^ (dashed lines) spectral func- 
tions in the Kondo case e d = -6A, T K /A = 8.07 x 1(T 4 (using T K = A 1 / 2 e~ 7red / 2A ), n f = 0.934. 
The + signs are for the spectral function above the threshold, Eq, and the circles are for the 
spectral function below the threshold. The arrow indicates the position of \ed\- 

FIG. 2. The exponents ccy^ deduced from the asymptotic power law behaviour of the auxiliary 
spectral functions as calculated within the NRG for different values of the occupation rif. The solid 
lines are the functions n/ — n"j/2 (o) and 1 — nj/2 (o). 
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